knitr::opts_knit$set(root.dir = '../../')
source('scripts/my_r.R')
library(feather)
## Warning: package 'feather' was built under R version 3.3.2
library(ggplot2)
library(reshape2)
library(bsselectR)
info <- c('Type1', 'newseq.matchgc.Noonan', 'motif_activation')
pos <- load_and_melt('pattern/Type1/newseq.matchgc.Noonan/motif_activation/pattern_positive.feather')
neg <- load_and_melt('pattern/Type1/newseq.matchgc.Noonan/motif_activation/pattern_negative.feather')
ggplot(pos$table.melted) + geom_raster(aes(x = variable, y = id, fill = value)) + scale_fill_gradient2() + ggtitle('Activation per sequence (Positive)') + labs(x = 'motif', y = 'sequence')
ggplot(neg$table.melted) + geom_raster(aes(x = variable, y = id, fill = value)) + scale_fill_gradient2() + ggtitle('Activation per sequence (Negative)') + labs(x = 'motif', y = 'sequence')
pvalues <- c()
motifs <- c()
for(i in colnames(pos$table)){
if(i != 'id'){
poss <- pos$table[[i]]
negs <- neg$table[[i]]
pvalues <- c(pvalues, wilcox.test(poss, negs)$p.value)
motifs <- c(motifs, i)
}
}
df <- data.frame(motif=motifs, p.value=pvalues)
##2 : Motif.291
##3 : Motif.275
##4 : Motif.100
##5 : Motif.59
##6 : Motif.78
##7 : Motif.241
##8 : Motif.29
##9 : Motif.168
##10 : Motif.193
##11 : Motif.73
##12 : Motif.46
##13 : Motif.106
##14 : Motif.0
##15 : Motif.154
##16 : Motif.143
##17 : Motif.222
##18 : Motif.278
##19 : Motif.133
##20 : Motif.14